addpath('./build')

load('a_opt');
load('x_opt');
x0 = x_opt;
a0 = a_opt;

access_A_matrix(a_opt);

p0 = [1 0];

numsteps = 1;

[xr pr] = reset_map(x0);
data = hybrid_flow(xr, pr, numsteps);

y = (a(1) * cos(a(2) * t) + a(3) * sin(a(2) * t)) / exp(a(4) * t) + a(5);